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Abstract 



A possibility for fundamental constants to vary in time is suggested by 
theories unifying gravity with other interactions. In this article we examine 
proposals to use optical transitions of Sr, Dy, Ybll and Yblll for the search 
of the time variation of the fine structure constant a. Frequencies of atomic 
transitions are calculated using relativistic Hartree-Fock method and configu- 
ration interaction technique. The effect of variation of a on the frequencies is 
studied by varying a in computer codes. Accuracy of measurements needed 
to improve current best limit on the time variation of a is discussed. 
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I. INTRODUCTION 



Theories unifying gravity with other interactions suggest that fundamental constants 
could vary in space-time (see, e.g. [1]). Recent evidence of variation of the fine structure 
constant a in quasar absorption spectra [2] elevated interest to the search of variation of a 
in laboratory experiments. Comparing frequencies of different atomic transitions over long 
period of time is a good way to do such search due to extremely high accuracy of mea- 
surements achieved for certain types of transitions. The best limit on local present-time 
variation of the fine structure constant published so far was obtained by comparing Hg + 
microwave atomic clock vs hydrogen maser [3]. Recently this limit was further improved 
by more than an order of magnitude in comparing cesium and rubidium atomic clocks [4]. 
There are also many proposals for the search of variation of a in atomic optical transitions, 
some of which were analyzed in our previous works (see [5] and references therein). In the 
present paper we analyze three new proposals involving strontium/calcium, dual beam [6], 
dysprosium atom [5,7] and ytterbium positive ions Yb + [8] and Yb 2+ [9]. We perform rela- 
tivistic many-body calculations to link variation of a with the variation of the frequencies of 
atomic transitions. Then we use this connection to find out what accuracy of measurements 
is needed to improve current best limit on time variation of the fine structure constant. 

In the proposal suggested by S. Bergeson strontium-calcium dual beam is to be used 
to compare the frequencies of the 1 S , — 3 P\ clock transitions in these atoms over a long 
period of time. Ca and Sr have similar electron structure. However, due to higher nuclear 
charge, relativistic effects are larger for strontium. If a is changing, corresponding change 
in frequency of the clock transition for Sr would go considerably faster than for Ca. Precise 
measurements might be able to indicate this or, at least, put strong constrain on possible 
variation of a. Calculations of the relativistic effects for Ca were done in our previous work 
[5]. In present paper we do similar calculations for Sr. 

Experiments with ytterbium positive ion have advantages of greater relativistic effects 
due to larger nuclear charge and the convenience of working with two different transitions 
of the same element. There are two transitions in Yb + involving metastable states for 
which comparison of frequencies is considered. One is quadrupole transition 4/ 14 6s 2 Si/ 2 — 
4f 14 5d 2 D 5 / 2 and another is octupole transition 4/ 14 6s 2 5*1/2— 4/ 13 6s 2 2 F 7 / 2 - The quadrupole 
transition is basically a s — d transition while the octupole one is a / — s transition. According 
to simple analytical formula presented in Ref. [5] relativistic energy shifts for s electrons, 
and electrons with high total momentum j (like d and / electrons) are large but have 
opposite sign. This means that we should expect that two metastable states of Yb + move 
in opposite directions if a is changing. This brings extra enhancement to the sensitivity of 
the measurements for Yb + to the variation of a. Our accurate calculations presented below 
support these considerations. 

The proposal for dysprosium is quite different from what was considered so far. Instead of 
comparing two very stable atomic clock frequencies the authors of this proposal [5] suggest 
to measure very small frequency of the transition between two almost degenerate states 
of opposite parity in dysprosium. The states are 4f 10 5d6s 3 [10]io E = 19797.96cm -1 and 
Af 9 5d 2 6s 9 Kiq E = 19797.96cm -1 . These states were used before for the search of parity 
non- conservation in Dy [11]. Small energy splitting and different electron structure of these 
two states lead to very strong enhancement of the sensitivity of the frequency of transition 
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between the states to variation of a. The enhancement (about eight orders of magnitude) 
seems to be strong enough to overcome the disadvantage of dealing with states which are 
not very narrow. 

In the present paper we calculate the values of relativistic energy shifts for Sr, Yb + and 
Dy and discuss what accuracy of measurements is needed to improve current best constrain 
on local time variation of the fine structure constant. 



II. CALCULATIONS OF ENERGIES 

We use relativistic Hartree-Fock (RHF) and configuration interaction (CI) methods to 
do the calculations. 

RHF Hamiltonian is used to generate a set of single-electron orbitals. We use a form of 
singe-electron wave function which explicitly depends on the fine structure constant a: 

*w,,. = - { . u Wt ) (d 

r \ iag ri {r)VL[n) flm J 
Then the RHF equation for ip(r) n has the following form (in atomic units) 

f' n (r) + ^fn(r) - [2 + a\e n - V HF )]g n (r) = 0, 

9n(r) ~ ^fn(r) + (e n - V HF )f n (r) = 0, 1 1 

here k = (— l)' +J+1//2 (j + 1/2), n is the principal quantum number and Vhf is the Hartree- 
Fock potential. The value of relativistic effects is studied by varying the value of a in (2). 
In particular, non-relativistic limit corresponds to a = 0. 

In the CI calculations we use approach similar to what was first developed in Ref. [12]. 
Electrons occupying open shells are considered as valence electrons and all other are core 
electrons. Correlations between valence electrons are treated within the CI method while 
correlations between valence and core electrons are included by means of the many-body 
perturbation theory. The effective CI Hamiltonian for valence electrons is presented by 

TV 2 

H CI = T.h, + T. e -- (3) 

i=l i<j V 

Here N is the number of valence electrons and hi is an effective singe-electron Hamiltonian 

Ze 2 

hi = cap + (/3 - l)mc 2 + V COTC + ^ (4) 

n 

Here V covc is the Hartree-Fock potential created by core electrons. It differs from Vhf in Eq. 
(2) by contribution of valence electrons. Si is so called "correlation potential" operator. It 
describes correlations between a particular valence electron and core electrons (see Ref. [12] 
for details). Note that in contrast with Ref. [12] we don't include in present work the S 2 
operator, which is a two-electron operator describing different type of correlations between 
valence and core electrons. Terms with S 2 can be considered as screening of Coulomb 
interaction between valence electrons by core electrons. These terms are less important 
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than those with Ei but much more time consuming in calculations. We either neglect them 
or simulate their effect by introducing screening factors. 

We are now going to discuss the specifics of the calculations for each atom/ion. Apart 
from the states of interest we also calculate energies of the other states of the same con- 
figurations to ensure that the accuracy is systematically good. We also calculate magnetic 
(^-factors to ensure correct identification of states. This is particularly important for dys- 
prosium. 

A. Strontium 

Strontium in its ground state is a closed-shell atom. It has two 5s-electrons on its 
outermost shell and we need to consider energy intervals between 1 S'o ground state and 
states of the 5s5p configuration where the 3 Pi metastable state is of most interest. The RHF 
calculations for Sr were done in V N approximation, for a closed-shell atom in its ground 
state. For the CI calculations we considered Sr as an atom with two valence electrons and 
followed the similar calculations for Ba [13]. Basis states for the CI+MBPT method were 
calculated using the B-spline technique [14] with 40 B-splines in a cavity of radius R = AOas- 
The same basis functions were used to calculate Ei and for the CI calculations. Thirteen 
lowest states above core in each of the S1/2, P1/2, P3/2, ^3/2 and (I5/2 waves were used to 
construct two-electron wave function for both 5s 2 and 5s5p configurations. Large number 
of basis functions is needed mostly for adequate description of the 5s5p configuration. This 
is because the V N approximation doesn't provide us with a good 5p single-electron state. 
Also, the 5s single-electron state in the 5s5p configuration is different from the 5s state in 
the 5s 2 configuration for which Hartree-Fock calculations were done. However, with thirteen 
states in each wave the saturation of the basis was clearly achieved and adding more states 
to the basis didn't change the energy. Two-electron basis states for the CI calculations were 
obtained by distributing valence electrons over 65 basis states (13 x 5) in all possible ways 
with a restriction of fixed parity and total momentum. 

The results are presented in Table I. As one can see the accuracy for the state of the 
interest 3 Pi is better than 1% while accuracy for other states is also good. 

B. Ytterbium 

The ground state of ytterbium positive ion is 4/ 14 6s 2 Si/ 2 and we need to consider 
transitions into the Af u bd 2 D 5 / 2 and 4/ 13 6s 2 2 F-j/2 states. Therefore it is convenient to do 
the RHF calculations in the V N ^ approximation, for the Yb 2+ ion with the 4/ 14 closed-shell 
configuration. The 6s, 5d and other basis states for the CI method are calculated then in the 
field of frozen closed-shell core of Yb 2+ . Then, in the CI calculations, we need to consider 
all 4/ electrons as valence ones since one of the transitions of the interest involves excitation 
from the 4/ subshell. So, the total number of valence electrons in present CI calculations 
is fifteen. This is very different from our previous calculations for Yb + [5] in which the 
4/ 13 6s 2 2 Ft/2 state was not considered and we were able to treat ytterbium ion as a system 
with one external electron above closed shells. 
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Our final set of single-electron states for the CI calculations consisted of 4/ 5 / 2 , 4/7/2, 
6S1/2, 5^3/2, 5^5/2 and few more s and / states above 4/ and 6s. Note that in contrast with 
Sr we don't need many basis functions here because all our single-electron wave functions 
correspond to the Yb + . This makes initial approximation to be very good and leads to fast 
convergence of the CI calculations with respect to the basis set used. 

We also don't include Si in calculations for Yb + . In a case of many valence electrons 
(fifteen for Yb + ) correlations are dominated by correlations between them which are taken 
into account accurately via the CI technique. Correlations between valence electrons and 
core electrons mostly manifest themself via screening of the Coulomb interaction between 
valence electrons. We take this effect into account semiempirically, by introducing screening 
factors fk- Namely, we multiply every Coulomb integral of the multipolarity A; by a numerical 
factor fk which is chosen to fit the energies. It turns out that good fit for Yb + is achieved 
with f 2 = 0.8 and fk — 1 for all other k. 

Many-electron basis states for the CI calculations were obtained by allowing all possible 
single and double excitations from the base configuration with the restriction of fixed parity 
and total momentum. 

Results for energies of Yb + are presented in Table I. The theoretical accuracy for energies 
as compared to the experiment is 2- 3% for the states of interest and is not worse than 5% 
for other states. 

C. Dysprosium 

Dysprosium atom is the most difficult for calculations because of its complicated electron 
structure. Ground state configuration of Dy is 4/ 10 6s 2 which means that there is no realistic 
RHF approximation which corresponds to a closed-shell system. We do the RHF calculations 
for Dy in the V N approximation with an open-shell version of the RHF method. Contribution 
of the 4/ electrons into the RHF potential is calculated as for a closed shell and then 
multiplied by a numerical factor to take into account its fractional occupancy. This factor is 
10/14 when interaction of the 4/ electrons with other core electrons is considered and 9/13 
when interaction of a 4/ electron with other 4/ electrons is considered. When convergence 
is achieved we have the 4/ and 6s basis states for the CI calculations. To calculate other 
states of valence electrons we remove one 6s electron, freeze all RHF orbitals, including 4/ 
and 6s and calculate the 6pi/ 2 , 6p 3 / 2 , 5d 3 / 2 , 5rf 5 / 2 and few more ci-states above 5d in the field 
of frozen RHF core. 

In the CI calculations states below 4/ are considered as core states and all other as 
valence states. Total number of valence electrons is therefore twelve. As for the case of 
Yb + we neglect Si and use screening factors as fitting parameters to improve agreement 
with experiment. It turns out that best fit for the 4/ 10 6s6p configuration is achieved with 
fx = 0.7 and fk — 1 for all other k. No fitting was used for other configurations. 

To calculate states of the 4/ 10 6s 2 , 4/ 10 6s6p and 4/ 10 6s5<i configurations we use the 
4/5/2, 4/7/2, 6si/ 2 , 6pi/2, 6^3/2, 5g?3/2 and 5d$/2 single-electron basis functions and all pos- 
sible configurations which can be obtained from these basis functions by exciting of one or 
two electrons from the base configuration. Same approach doesn't work for the Af 9 5d 2 Qs 
configuration because of huge number of many-electron basis states generated this way and 
as a consequence, the CI matrix is of so large size that it could not be handled by our 
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computers. On the other hand test calculations with pairs of configurations showed that 
mixing of our state of interest with other configurations is small and can be neglected. We 
do need however to include mixing with the 4f 9 5d6d6s, 4f 9 5d7d6s and 4f 9 6d 2 6s configura- 
tions. This is because our basis 5d state corresponds rather to the Af 10 5dQs configuration 
and extra d-states are needed to correct it. 

The result are presented in Table I. Note that they are considerably better than in our 
previous calculations [15]. This is because of better basis and more complete CI treatment. 



III. FREQUENCY SHIFT: RESULTS AND DISCUSSIONS 

In the vicinity of the physical value of the fine structure constant (a = a ) frequency 
(cu) of an atomic transition can be presented in a form 

00 = uj + qx, (5) 

where x = (a 2 / a o) — 1, cu is the experimental value of the frequency and q is a coefficient 
which determines the frequency dependence on the variation of a. To find the values of q for 
different atomic transitions we repeat all calculations for a = a = a and a = \J~^cn . 

Then 

q = A{u + -uJ), (6) 
Where u + is the value of u for a = J f a , an d W- is the value of uo fora = J ^a . Calculations 



for a = a are done to compare the theory with experiment and to check whether frequencies 
are linear functions of a 2 . The results for coefficients q are presented in Table II. Note that 
we have included in the Table the results of our old calculations for Yb + . These calculations 
were done in a very different way, assuming that Yb + is an atom with one external electron 
above closed shells. Comparison of the results obtained by different methods gives estimate 
of the accuracy of calculations. 

Search for the time variation of the fine structure constant can be conducted by compar- 
ing two frequencies of atomic transitions over long period of time. The measured value can 
be presented as [5,3] 

AW = (±-±). (7) 



Using Eq. (5) one can reduce Eq. (7) to 



A( t ) (A). (8) 

Current best laboratory limit on the time variation of a is a/a < 10~ 15 yr^ 1 [4]. 

In the first experiment considered in this paper a dual calcium-strontium beam is to be 
used to compare the frequencies of the 1 S'o — 3 Pi transitions in both atoms. Substituting 
Ui = 15210cm -1 , q 1 = 230cm~ 1 for Ca [5], u 2 = 14504cm" 1 , q 2 = 667cm" 1 for Sr (Tables 
I, II) and a/a — 10" 15 yr" 1 we get 
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A(t)(Sr-Ca) = 6.2 x 10 _17 yr _1 - 



(9) 



Note that the width of 3 Pi state in Sr may be a problem in this case. 

In the case of Yb + frequencies of the 2 Si/ 2 — 2 -D5/2 and 2 Si/ 2 — 2 P7/2 are to be compared. 
Substituting the numbers we get 

A(t)(Yb + ) = 6.1 x Kr 15 yr _1 - (10) 

Note two orders of magnitude improvement in the magnitude in comparison with the Sr-Ca 
dual beam experiment. 

We have also calculated (/-coefficient for 4/ 14 1 S'o - 4f 13 5d 3 Po (w = 45276 cm' 1 ) tran- 
sition from Yblll ground state. This was motivated by the proposed measurements [9] of 
a-variation using comparison of 1 S - 3 P transition frequencies in In + , Tl + and odd isotope 
of Yb ++ . The different signs and magnitudes of relativistic corrections in In + (q= 4414 
cm -1 ), Tl + (g=19745 cm -1 ) and Yb ++ (q= -27800 cm -1 ) provide an excellent control of 
systematic errors since systematic errors are not correlated with signs and magnitudes of 
the frequency shifts qx, where x = (a 2 /al) — 1. The same idea (combination of anchors, 
positive shifters and negative shifters) has been used to control systematic errors in Ref . [2] . 

In our view, a very interesting possibility is that for dysprosium. Instead of comparing 
frequencies of different transitions one should measure the energy difference between two 
very close states of opposite parity. The corresponding g-coemcient is q = 6008 + 23708 = 
29716cm -1 (see Table II). The frequency of this transition ranges from few MHz to few GHz 
depending on isotopes and hfs components used. If we take, e.g. uo = 3.1MHz [10] we get 

A(t)(Dy) = 5.7xl0 8 (-). (11) 

a 

This is an eight orders of magnitude enhancement in the relative value of the effect compared 
to atomic clock transitions! Substituting a/a = 10 _15 yr _1 we get 

A(t)(Dy) = 5.7 x 10 - V -1 - ( 12 ) 

This means that to improve current best limit on local time variation of a the frequency of 
this transition in Dy should be measured to the accuracy of about 10~ 7 over about a year 
time interval. This seems to be feasible [7]. 
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TABLES 



TABLE I. Energies and g- factors 



Atom/ion 


State 


Energy, 1/cm g 






theory experiment theory experiment 3, non-relativistic 



5s 2 


'So 














5s5p 


3 Po 


14171 


14318 








5s5p 


3 Pi 


14384 


14504 








5s5p 


3 P2 


14832 


14899 








5s5p 


'Pi 


22829 


21698 








4/ 10 6s 2 


5 h 








1.243 


1.242 


1.25 


4/ 10 6s 2 


5 /7 


4123 


4134 


1.175 


1.173 


1.179 


4/ 10 6s 2 


5 /6 


7147 


7051 


1.072 


1.072 


1.071 


4/ 10 6s 2 


"h 


9428 


9212 


0.907 


0.911 


0.9 


4/ 10 6s 2 


5 /4 


11199 


10925 


0.614 


0.618 


0.6 


4/ 10 5ti6s 


3 [8] 9 


18605 


17515 


1.319 


1.316 




Af 10 5d6s 


3 [9]io 


18615 


18463 


1.291 


1.282 




4f 10 5d6s 


3 [10]n 


19811 


19349 


1.268 


1.27 




Af 10 5d6s 


3 [10]io 


20133 


19798 


1.208 


1.21 




4/ 9 5d 2 6s 


9 K 12 


23345 


22541 


1.327 


1.333 


1.333 


4/ 9 5d 2 6s 


9 ^n 


20513 


20448 


1.352 


1.354 


1.303 


4/ 9 5ci 2 6s 




19623 


19798 


1.372 


1.367 


1.264 


4f 9 5d 2 6s 


% 


19434 


19558 


1.390 


1.39 


1.377 


4/ 9 5d 2 6s 


9 G 8 


18379 


18473 


1.461 


1.46 


1.5 


4f 9 5d 2 6s 


9 G 7 


18662 


18529 


1.492 


1.467 


1.5 


4f 9 5d 2 6s 


7 F 6 


19714 


19304 


1.527 


1.54 


1.5 


4/ 9 5d 2 6s 


7 G 5 


21697 


20892 


1.510 


1.32 


1.37 


4/ 9 5d 2 6s 


9r< 


23748 


22697 


1.492 


1.487 


1.5 


4/ 14 6s 


2 $l/2 








2.000 


1.998 


2 


4/ 14 5d 


2 F>3/2 


22888 


22961 


0.800 


0.800 


0.8 


4f 14 5d 


2 D 5/2 


23549 


24333 


1.200 


1.202 


1.2 


4/ 13 6s 2 




31820 


31568 


0.857 


0.862 


0.857 


4/ 13 6s 2 


2 F 7/2 


21819 


21419 


1.143 


1.145 


1.143 


4f 1A 6p 


Pl/2 


26000 


27062 


0.667 


0.667 


0.667 


4f 14 6p 


2 P 3/ 2 


29005 


30392 


1.333 


1.333 


1.333 



a References [16,17] 
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TABLE II. Relativistic energy shift q (cm 1 ). 



Atom/ion 


State 




This work 


[51 


Sr 


5s5p 


3 i>! 


667 






5s5p 


^1 


1058 




Dy 


Af l0 5d6s 


3 [10]io 


6008 






4/ 9 5d 2 6s 


9 K 10 


-23708 




Ybll 


4/ 14 5d 


2 D 3/2 


10118 


12582 




4/ 14 5d 


2 D 5/2 


10397 


11438 




4/ i3 6s 2 


F 7/2 


-56737 




Yblll 


4/ 13 5d 


3 Po 


-27800 





